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We begin by introducing a new procedure for construction of the exact solutions to Cauchy 
problem of the real-valued (hyperbolic) Novikov-Veselov equation. The procedure shown therein 

utilizes the well-known Airy function Ai(^) which in turn serves as a solution to the ordinary 
,2 

differential equation In the second part of the article we show that the aforementioned 

procedure can also work for the n-th order generalizations of the Novikov-Veselov equation, provided 
that one replaces the Airy function with the appropriate solution of the ordinary differential equation 


PACS numbers: 


I. INTRODUCTION 

The Novikov-Veselov equation (NV) first came into being almost 30 years ago, in 1984, when Sergey Novikov and 
Aleksander Veselov have introduced it as a two-dimensional version of the Korteweg-de Vries equation [NV]. Since 
then, a huge and ever-growing body of works related to the study of NV equations has been established (see, for 
example, [GN], [LMSS], [P] and also [CMMPSS] for a rather extensive review of a recent literature on the subject). 
One of the most prominent aspects of this interest was an almost universal adoption of an inverse scattering method 
as a primary tool for conducting the research and finding the exact solutions of NV. However, in this article we wish 
to discuss an alternative method of solving the Cauchy problem for NV. This method, albeit simple in principle, 
appears to be deep enough to be applicable to a rather broad class of equations, NV and the 2-dimensional nonlinear 
Schrodinger equation being its two prominent members. 


II. THE MOUTARD TRANSFORMATION 


Let us start by introducing the hyperbolic NV equation: 


Ut — tlxxx ^yyy “b -p (bu')y^ 

ttx — tly, Uy = bx- 

This system allows for a Lax pair of the following type: 


( 1 ) 


Aixy -P = 0 

^XXX “P ^ yyy + 3{a'4/x + b'iiy). 


( 2 ) 


If one knows two linearly independent solutions 4'i(a;,j/,t) and 4 / 2 ( 2 ^, y, t) for (2), then one can utilize the famous 
Moutard transformation to construct a new function 'k[l](a;, y, t) that will serve as a solution to the same equation 
(2) albeit with a new potential u[l](x,y,t). The new potential will then satisfy the relation 

u[l] = u + 2dxdy\D.A!i. (3) 

Let us assume that u = a = b = 0. Then the entire system (2) simplifies to 


0 

II 

(4) 

= 'ilxxx + '^yyy 

(5) 
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The equation (4) can be resolved by separating the variables. The resulting solution will be of a form: 

= A[x,t) + B{y,t), (6) 

where A^B are two arbitrary functions that are continuously differentiable by x and y, correspondingly. Substituting 
(6) into (3) yields a following post-Moutard form of function u[l]{x,yA)'- 


, , dxA ■ dyB 

"W = 


( 7 ) 


As follows from (7), our next goal should lie in ascertaining the exact forms of the functions A{x, t) and B{y, t). This 
task can be accomplished by looking at the equation (5) which we have ignored so far. We will rewrite it as a standard 
Cauchy problem by introducing the initial conditions for A(t, x), B(t^ y) 


and rewriting the (5) as a system 


A{0, x) = (j){x), B{0, y) = $(?/). 

At = Axxx + 

Bt = Byyy “ T (t) , 


( 8 ) 


(9) 


where T = Tit) is an arbitrary time-dependent function. The apparently symmetric nature of (9) allows us to restrict 
our attention on just one of the equations therein, namely - the first one. 

We begin by introducing the Fourier transform A(j),t) of the function A(xA)'- 


A(p, t) = 


1 




A{x,t)e-^P^dx. 


This transformation is handy because of the identity 

OO 

— OO 

which, after being substituted into (9), yields the equation 


e^P^dp, 


dA 

'm 


ip^A - T e*P^dp = 0. 


( 10 ) 


( 11 ) 


The equation (11) must be satisfied for all x and p, and therefore leads to: 

(12) is a nonhomogeneous linear O.D.E. of first order. Its general solution is 


Z 

A{pA) = C{p)e-^p"^ + j T(r)e-*P'(*-^)dT, 


( 12 ) 


(13) 


where C{p) is a function, determinable from the initial conditions (8). Using the inverse Fourier transform (10) we 
come to the following conclusion: 


OO / t > 

A{x,t) = ^ J I J T{T)e^p'-dT + C[p) 


ipx — ip^t 


dp. 


(14) 


-OO \0 


According to (8), 


0(x) = 




C{p)e^P^dp, 


— OO 
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so the unknown C{p) is an inverse Fourier transform of the initial condition i.e.: 

CO 

C{p) = ^ / <i^{x)e-^^^dx. 

— CO 

and we subsequently end up with the following general formula for the function 

t CO CO CO 

A{x,t) = -^ J T{T)dT j ^ J J (15) 

0 —oo —oo —oo 

The (15) can be further simplihed by pointing out the similarity between the integrals with respect to variable p 
and the Airy function Ai(^). The Airy function is a particular solution of the eponymous Airy equation: 


d'^z 


2 


that has a following integral representation: 


OO 




dt. 


(16) 


(17) 


Using this fact together with the apparent identity: 


1 


^dp = 


1 


Ai 


. 1 - 


— CO 

together with the equation (15) and the similar one written for B{y,t) finally yields: 


I 

A{x,t) = J 


B{y,t) = - 


T{t) 


Ai 


\/3(t - 1) 


dr - 


1 


r(r) 


Ai 


^3(r - t) I ^3(r - t) 


dr 


1 




</'(OAi df 


^{p)Ai ( ^ ^ ) dr]. 




(18) 


So, we end up with both the solution = A + i? of the Lax pair (4), (5), and, as a courtesy of Moutard transform 
(3), with a solution m[ 1] of the NV equation (1) as well. In other words, to hnd a non-zero solution of the NV 
equation, it will suffice to start with m = 0, impose the boundary conditions (8) on the Lax pair (4), (5), use (18) 
to find its solution and conclude the calculations by finding a function m[1] via the Moutard transformation (3). As 
straightforward as it is, there is one question we should ask: what would happen should we try to invert the process 
and instead start out with he boundary conditions for the NV equation itself? 


III. THE CAUCHY PROBLEM FOR THE NOVIKOV-VESELOV EQUATION 

In the previous chapter we have shown that there shall exist a solution u[l](a;, y, t) to the NV equation, whose exact 
form can be derived via the Moutard transformation (7) from the solutions of the system (4, 5), provided we are given 
the initial conditions (8). But what would happen if the exact forms of the functions (fix) and ^{y) are unknown and 
we are instead given the initial conditions for the NV equation itself, and would it still be possible to find the required 
m[ 1]? In other words, is it possible to find an analytic solution to the Cauchy problem for the NV equation provided 
we only know that the solution has a general structure (7)? In short, the answer is “yes”. 

Let us start by introducing the set of conditions and boundary conditions for the NV equation: 

u[l]{x,y,0) = uo{x,y) 
uo{x,0) = Ai{x) 
uo(0,y) = Biiy) 

Ai(0) =Bi(0) = C, 


( 19 ) 
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where C G M is some constant that is given to us together with the boundary conditions Ai and -Bi. Since we know 
that m[ 1] satisfies the Moutard transformation, we also know that: 


uo{x,y) = -2 


(j)'{x) ■ 4>(y) 


( 20 ) 


where cj) and $ are defined as in Sec II, and ' and • denote the partial derivatives with respect to x and y variables 
correspondingly. From (20) and (19) it immediately follows that 


Ai{x) = —2 


B,{y) = -2 


(j)'{x) ■ $(0) 

(<)>(x) + $(0))2 

■ $(z/) 

( 0 (O) + $(y ))2 


p g <^'(0) ■ m 


( 21 ) 


The first two differential equations in (21) can be easily integrated; for example, the first one after the integration 
with respect to the variable x yields 


^(0) 

cj){x) + $(0) 


i I A,{Od^ + 

0 


^(0) 

(^(0) + $(0)’ 


which leads us to the following conclusion: 


(j){x) 

0'(x) 


2^0 

f Ai{^)d^ + 

0 


$0 


-24>o^i(a^) 




2 ’ 


( 22 ) 


where we have introduced the notation: (/)(0) = (j)o, <I)(0) = <I)o, ^'(0) = and 4)(0) = 4>o. In a similar fashion, the 
boundary condition $(?/) and its derivative will satisfy: 


^(y) 

^(y) 


2</>(, 


jB^iOdC 


201, 

<?^0 + ' J ’0 


— <)>o 


-2^oBi{y) 


/Bi(C)dC + ^^^ 


2 ■ 


<pQ-\-^Q 


(23) 


The system (22, 23) depends on four constants: (po, <I>o, pQ and 4)0. Three of them can be chosen arbitrarily, whereas 
the fourth one would have to satisfy the equation (21), namely: 

g 4>0 ■ ^0 ^ 

{po + <l)o)^ 

Curiously, this choice does not affect the Cauchy problem of the NV equation in the least, for it can be shown by 
direct substitution into (20) that: 


uo{x,y) 


4CAi{x)Bi [y) _ 

' X y \ 2 ’ 

ffA,(0B,(0d(:d^ + 2C] 

.0 0 J 


(24) 


i.e. the initial condition uo{x,y) depends only on the known initial boundary conditions Ai{x)^ Bi{y) and C. 
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We are now ready to answer the question posed in the beginning of this section: provided we know the initial 
conditions (24), how do we solve the corresponding Cauchy problem of the hyperbolic real-valued Novikov-Veselov 
equation? The answer lies in repeating the Moutard transformation process we described in Sec. II! Indeed, since 
the unknown functions (f>{x) and 4>(y) satisfy the relations (22) and (23), all we really have to do is substitute them 
into the system (18), derive A{x,t) and B{y,t), and substitute them in equation (7) to find out the sought after 
u[V\{x,yA), which will conclude the problem. 

Lets summarize everything we have said so far. In order to find an exact solution u(x,y,t) to the hyperbolic 
real-valued Novikov-Veselov equation 


Ut 


^yyyi 


(25) 


with the given initial boundary conditions 

u{x, 0 , 0 ) = Ai{x), 


u( 0 ,y, 0 ) = i?i(y), u( 0 , 0 , 0 ) = C, 


that correspond to the initial condition 


uo{x,y) 


ACAi{x)Bi{y) 


JJA^i0B,{0dCd^ + 2C 


2 ’ 


one shall: 

1. Choose a differentiable function T(t) and four numbers a, /3, 7 and S that satisfy the condition, 

7 • S 


-2 


{a + l3y 


= c. 


2. Find two support function 4>{x) and 4>(y) via the formulas 

- 


$(y) = 


f Ai(^)d^ + ^jTa 
0 

27 

fBi(C)d(+^ 


— a. 


3. Substitute and ^(y) into the equations 


t 

A(x,t) = J 


T{r) 

^3(t — t) I -y3(T — t) j 


Ai 


dr 




B{y,t) = - f 


V ^3(T-t) 


dr - 




$(77)Ai 


■n-y 


drj. 


(26) 


4. Substitute the new functions A{x,t) and B{y,t) into the equation 


u = —2 


dxA ■ dyB 

{A + BY ■ 


(27) 


The resulting function u{x, y, t) will be a proper solution of the Cauchy problem since by construction it will satisfy 
both the NV equation (25), and the initial conditions m(x, y, 0) = 'Uo(x, y). We would like to emphasize here that this 
procedure does not involve anything more complicated than partial differentiation and integration and can therefore 
be used for both the analytic study of the properties of the solutions of NV equation and the corresponding numerical 
calculations. 
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IV. GENERALIZATION OF THE METHOD: THE HIGHER ORDER EQUATIONS 

Let us now say a few words about the more general problem. Suppose we have the following operator-type Lax 
pair: 


dxdy^ + wdr = 0 


where n G in some non-zero natural number. This system will correspond to a family of Lax equations, with the 
special case n = 3 corresponding to the hyperbolic NV equation. It will still allow for the Moutard transformation, 
and therefore the crux of our discussion would still be applicable for arbitrary n. However, one thing that must 
change is the exact form of the equations for A{x,t) and The formula (26) is no longer applicable for the 

general case and should be properly replaced. In order to find out the suitable replacement, we shall separately 
consider two alternative cases: when n is odd and when n is even. 

Case 1: Odd n. Let n = 2m + 1, where m > 0. Following our previous discussion, let us consider the special case 
M = 0. Then the system (28) turns into 


dxdy'i' = 0 

Since the first equation requires that d' = A(x, t) + B{y, t), the system subsequently splits into the following equations: 

dtA2m+l = dtB2m+l = B2m+1, 

where for simplicity we have omitted the arbitrary function T{t). Using the Fourier transformation 

OO 

^2m-t-l(P;t) — y/2 — j‘ ^ dx, 

— OO 


we end up with the differential equation 

dtA2m+l = {ip)^^+^A2m+l = ^(- 1 ) • 

Solving it and returning back to A(x, t) as described in Sec. II yields 

OO 

I 


A2m+i{xA) = ^ I I 


(29) 


As we know, in the special case m = 1 (i.e. n = 3) the inner integral in (32) can be rewritten in terms of the Airy 
function 


= wJ‘' 

— OO 

which serves as a solution to the Airy equation 


*3 +u) ^ J cos dt, 


d^z 


and is easily derived using either Fourier or Laplace transformation [6]. 

Similarly, it is easy to show that one of a solutions to a more general equation 


d^2m 


= 
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will be a higher-order generalization of the Airy function [7]: 


Ai 2 m+l(C) — 




l( + 2 Tn.+ l'\ / 2 

2 m+l * J_ 


cos ( - I f 2 m+l 

2 m + 1 


jdt, 


(30) 


which means that the required functions A 2 m+i B 2 m+i can be derived from the initial conditions (j^ix) and $(y) 
by the following formulas: 


A 2 m+l(^i^) — 


+^( 2 m + l)t 


HO Ai2m+1 


^-x 


— OO ' 

OO / 

“ -y(2m + l)t / ( “ 


+■{/ (2m + l)t 


(31) 


C- y 


V( 2 m+l)t 


Case 2: Even n. Let n = 2m, where m > 0. This time let us utilize not a Fourier but a Laplace transform: 

OO 

A 2 m{x,t) = j A 2 m{p,t)e^'°dx. 

— CO 

the equation for A(p, t) is 


dA 


2m _ 2m 


dt 


— P -^2m 5 


so the required function Aix^ t) will satisfy the equation 

OO OO 

A2m{xA) = J d^ J dp 


(32) 


It is not difficult to show that the Laplace transformation method applied to the ordinary differential equation 

d 2 m-i^ ^ 

C2m-1 ~ 


dH 


will yield a following solution 


^ 2 

Ai 2 m (0 = J exp dt, 


(33) 


and so the even case produces the formulas that are quite similar to the old ones, namely: 


A2m{x,t) = 

B2m{x-i 0 — 


1 


^^/—2mt 

1 

^y—2mt 


di HO Ai2m 


dC ^>(0 Ai2m 


'V- 2 mt 

y-C 
‘^\/—2mt J 


(34) 


Note the appearance of a negative sign under the root in (34), which serves as a indication of an ill-posedness of our 
problem for t > 0. 
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